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General asymptotic approach to the stability problem of multi-parameter solitons in Hamiltonian 
systems idEn/dz — SH/5E^ has been developed. It has been shown that asymptotic study of the 
soliton stability can be reduced to the calculation of a certain sequence of the determinants, where 
the famous determinant of the matrix consisting from the derivatives of the system invariants with 
respect to the soliton parameters is just the first in the series. The presented approach gives first 
analytical criterion for the oscillatory instability and also predicts novel stationary instability. Higher 
order approximations allow to calculate corresponding eigenvalues with arbitrary accuracy. 
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T-H ■ 
bX)' I. INTRODUCTION 

'^~< ' Solitary waves ('solitons') can appear when an initial excitation applied to a medium is strong enough to cause 

^\ : nonlinear response. Formally solitons are solutions of some nonlinear partial differential equations and their dynamics 
T-H generally is a complex phenomenon, which can be described exactly only in the very special integrable situations 

y. The problems of soliton stability and instability induced dynamics in nonintegrable Hamiltonian models have 
\^ ' paramount importance for understanding of a wide range of physical phenomena covering such fields as propagation of 
J^ I electromagnetic, water and plasma waves, condensed matter physics and classical field theory |g-p. Several analytical 
r-^ ■ approaches to the stability problem are known. For instance, in the nearly integrable situations the perturbation theory 
f~^ ■ based on the inverse scattering transform can be used R,0| . Far from the integrable limit variety of methods can be 
\^ . applied. Among them asymptotic stability theory |0,|| , method of adiabatically varying soliton parameters ||l|,p|-[l3[ , 
^^ ' Lyapunov Q and Evans Q methods. 
^^ , Generally, stability of a solitary wave in a Hamiltonian model can be lost due to bifurcations involving appearance 

of a positive eigenvalue {stationary instability) in the soliton spectrum or a pair of complex conjugate eigenvalues 
T^ ] with positive real parts {oscillatory instability) |14| . Both types of these instabilities have been extensively studied in 
r^ ' the different solitonic contexts proving their ubiquitousness and fundamental importance, see Refs. |q-Mjlq,nq] and 
I ] Refs. p, pT] , p^ , p7[ , respectively, for the stationary and oscillatory instabilities. In the most of the known cases the loss 

of stability is associated with the collisions of the purely imaginary eigenvalues corresponding to the so called internal 

modes 20 1 of the soliton spectrum (see [p| for interesting exceptions). 
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Applying the above mentioned methods it was shown that in many cases a threshold of stationary instability of 

^ • multi-parameter solitons is given by the zero of the determinant of the Jacobi matrix Jij — d^^-Qi, where kj are the 

k>( , soliton parameters and Qi are the associated motion integrals [p|- [ll|jl^ ] . The condition det{Jij) = is, in fact, the 

j_] ■ compatibility condition of the problem arising in the leading (zero) order of the asymptotic solution of the eigenvalue 

Cd ' problem governing stability of the soliton |8| |l(]|Jl5|| . To find expressions for the eigenvalues it is necessary to proceed 

further and solve problems arising in higher (at least first) orders. Up to now this was done only for the specific 

class of model equations having single parameter soliton families [[7|J^] . For stationary bifurcations of two-parameter 

solitons adiabatic method has been applied in Refs. | pO| , pTp^ ]. Linear approximation of this method actually gives 

an expression for eigenvalues, for more details see Section IV. However, all known developments of this method fail 

to give criterion indicating transition to the oscillatory instability, i.e. instability with complex eigenvalues. It is also 

difficult to extend this method beyond its first order because of the rather involved calculations. 

The purpose of this work is to formulate a general asymptotic approach to stability of multi-parameter solitons in 
Hamiltonian models, to show how it can be used to find expressions for the instability growth rates with arbitrary 
accuracy and to formulate criterion for the oscillatory instability of solitons. 
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II. MODEL EQUATIONS AND SYMMETRIES 



We will consider Hamiltonian equations in the form 

.dEn _ 6H 

*~a7 ~ Ie* 



,n=l,2...N, (1) 



which describes a wide range of physical phenomena related with self-action and interaction of slowly varying wave 
envelopes in a variety of nonlinear media [|l&-|l2,n3-n^ , for general review of the Hamiltonian formalism see p8| . 
Here E„ are complex fields, z the propagation direction of the interacting waves, x the coordinate characterizing 
dispersion or diffraction, H = H{dxEn, En, dxE^, E^) is the Hamiltonian and * means complex conjugation. We will 
assume that H is invariant with respect to a set of (i — 1) phase transformations: 

En ^ En eMnnM, ^ = 1, 2, . . . (i ~ 1) , (2) 

(j)i are arbitrary real phases and jni are some constants. Because H does not depend on x explicitly, Eqs. d3) are also 
invariant with respect to arbitrary translations along x: 

En{x,z) ^ En{x ~Xo,z). (3) 

Symmetry properties (y), (^ together with Hamiltonian nature of our problem imply presence of L conserved 
quantities, see, e.g., p9[ , which are the {L — 1) energy invariants 

N 

Qi= dxJ2jni\En\\ / = 1,2,...(L-1), (4) 

•' n=l 

and the momentum 

N 

Ql^I^ dxY, {EldxEn - EndxEl). (5) 

* -^ n=l 

Another important consequence of the invariances (g), (^ is that a certain class of solutions of Eqs. (Q) can be 
sought in a form when xq and 0; are linear functions of z, i.e. Xq = klz and (j)i = kiz, then 

En{x, z) = a„(x - klz) exp(i ^ J^iiniz), (6) 

1=1 

where {^jj^J^ and kl are real parameters characterizing, respectively, phase velocities of the interacting waves and 
the soliton group velocity. Functions an{T) obey a system of the ordinary differential equations 

{iKLdr + a„)a„ = — —^, (7) 

where Ha = H{dTan, an, c^rO* , a*), t = x — klz and a„ = X]i=i imm- We assume now that in a certain domain of 
the parameter space (ki, k.^,- ■ ■ hl) Eqs. (0) have a family of the solitary solutions such that |a„| — *■ for r -^ ±oo. 

III. ASYMPTOTIC STABILITY ANALYSIS 

To study stability of the solitons we seek solutions of Eqs. ([|) in the form 

L-l 

En = {an{T) + £„(t, z)) exp(i ^ JnlKlz), (8) 

1=1 

where e„(T, z) are small complex perturbations. Linearizing Eqs. (n|) and assuming that £n(T, z) = £,n{T)e^^ , ejJjlT, z) — 
^„+Ar(r)e^^ we get the following nonself adjoint eigenvalue problem (EVP) 



where ^ = (^i, . . . ^^r, ^at+i . . . (,2n)'^ , and R, S are N x N matrix operators with elements given by 

r / , . c, \ , ^ ^a ^ S Ha 

Snl = Onl[an + JKid^j + , r„; 



(5a* (5a/ ' (5a* (5a;* ' 

here J„i is the Kroneker symbol. Note, that the operator S" is a selfadjoint one, i.e. S ^ S\ and i? is a symmetric 
operator, i.e. R — R^ . 

To solve EVP (|9|) we apply the asymptotic approach, which relies on expansion of the unknown eigenvector ^ into 
an asymptotic series near either neutral eigenmodes Qfll, i-e. zero-eigenvalue modes, of the operator £, or modes 
of continuum pG] , or both of them |jl2] . The neutral modes can be generated by infinitesimal variations of the free 
parameters of the soliton and thus always be presented as explicit functions of the soliton solution. At the same time 
continuum eigenmodes are explicitly known in the very rare, normally in integrable, situations jl2J2fl]. This is an 
important fact which makes the asymptotic expansion near the neutral modes by the very practical tool of the stability 
theory. However, as any approximate method, it has a certain limitation. Namely, it describes only eigenvalues A 
corresponding to a specific class of the perturbations which in the zero approximation can be expressed as a linear 
superposition of the neutral eigenmodes. Thus, generally speaking, on the basis of this approach one can get only 
sufficient conditions for soliton instability or, in other words, necessary conditions for soliton stability. Therefore 
presence of other instabilities which can be captured only numerically can always be expected |16|] . 

By infinitesimal variation of (f)i and xo it can be shown that 

ui = {jiiai,. ..jNiaN,-liia*i,.- ■ -Inio^nY , ul = Tp, 

OT 

a = (ai, . . . ajv, a*, . . . a^)^, I ^ 1,...{L-1) 

are neutral modes of C, i.e. Cui — {I = 1, . . . L). C also has L associated vectors Ui = da/dni such that 
CUi = -Ul, I ^ 1,. . .L. 

It is straightforward to see that any solution of EVP (g) must obey L solvability conditions 

{wi\XO=0, l = l,...L, (10) 

where {y\z) = X]i=i I ^^xy^Zi and wi are the neutral modes of the operator C\ C^wi — 0, 

wi = ['^iiai,. . ."fNia.NniiO'i^- ■ -IniOn) , WL = i-^, 
b = (-ai, ... - aN, al,. . . a^) , Z = 1, . . . (L - 1). 

Associated vectors of C'^ are Wi = db/dm and they obey OWi = —wi, I = 1, . . . L. 

Close to instability threshold it is naturally to assume that |A| ~ e ^ 1. As it was already discussed above we 
will consider a special class of the perturbations which in the leading approximation can be presented as a linear 
combination of the neutral modes. Therefore we seek an asymptotic solution of EVP (^ in the following form 

oc L 

m=0 1 = 1 

where constants Ci and vector-functions ^m>o have to be defined. Here and below I — 1,2, ... L. Substitution (O) 
into EVP M) gives a recurrent system of equations for ^m 



im>0 



a^_, 



e 



Co. (12) 



Substituting (|ll|), (O) into conditions (hOl) one will find the homogeneous system of the L linear algebraic equations 

oo L 



m=0 1=1 



for L unknown constants Ci. System (|3|) has a nontrivial solution providing that the corresponding determinant is 
equal to zero. This determinant is an infinite-order polynomial with respect to A^, which, in fact, is the asymptotic 
expansion of Evans function B . Zeros of this polynomial define the spectrum of the solitary wave linked with the 
chosen class of the perturbations. Thus the equation specifying eigenvalues A is 



A^^V(-A^F^,=0, (14) 






where Dj are the real constants. Eq. (14) always has zero root of the 2i-order. It indicates that each of the zero 
eigenvalues corresponding to the neutral modes ui is doubly degenerate one. This degeneracy originates from the 
presence of the associated vectors Ui . 

To write the explicit expressions for Dj it will be convenient to introduce vectors M™ = (A^J™ . . . Al;™ ), where, 

mP = {wi\C-^"-^U^,), m = 0, 1, . . . cx). 
Now each Dj can be presented as 

D,= Y. 2?(A^l"^\...A^f^^), (15) 

mi+...mL=j 

where I?(A^j^™ , . . .M™^) is the determinant of the L x L matrix consisting of the rows M™ and the sum is 

taken over all such combinations of (mi, . . . ttil) that X]/=i ^^i ~ J- -^ii' *^^^ ^^ readily expressed via derivatives of 
the conserved quantities with respect to the soliton parameters: 

and practical calculation of Al T for to > can be simplified: ^A\? = — {Wi\jC'^^~^"^^Ui') . Note, that in most of the 
cases solitary solution itself can be found only numerically using any of the well established methods for solving the 
nonlinear ode's. Recurrent calculations of CS^^'^'^'lJi can be readily reduced to the numerically even simpler problem 
of solving of the linear inhomogeneous ode's. 

Because |A| was assumed to be small, Eq. (141) has an asymptotic character. Therefore to make it work some 
additional assumptions must be made about orders of Dj . If these assumptions are satisfied then Eq. ( p^ describes 
correctly the soliton spectrum and predicts bifurcations of the soliton. The corresponding eigenvalues can be found 



using Eq. (14) with any degree of accuracy. For example, let us assume that Dq ~ e^ and -Dj>o ~ 0{1). Then, 
presenting /ras 

oo 

A^ = ^'E0'G■-e2^ (17) 

in the first order Eq. (|4|) gives a linear equation for Co, 

Do - e\oDi = 0, (18) 

which indicates a threshold of the stationary bifurcation at Dq = 0. This is precisely the condition det{Jij) — 
discussed in the introduction. Continuing to the next order one obtains 

A) / _ DoD_2 



X' = ^[l-^ + 0{e-)]. (19) 



If Dl ^ e^ then the asymptotic expressio n (|l9| ) fails. To have a balanced equation for C,q, we must now assume that 
Dq ^ e^ . However, in such a case the Eq. QSlfor Co changes from linear to quadratic: 

Do~e^Q^Di+e\lD2^Q. (20) 

Eq. ( po| ) gives two threshold conditions D^ ~ and Df = ADqD2, see Fig. n^. The latter condition indicates onset of 
the oscillatory instability for 



Dl < 4DqD2 



(21) 



Thus we have formulated analytic criterion for the oscillatory instability. It is also clear that the point Dqi = is a 
source for the novel stationary instability, see rightmost region Df > 4Z?o-D2, Di > in Fig. 11], where an eigenvalue 
which is positive throughout this region can not be predicted by Eq. (|l9|). 

It follows by recurrence that if £^,'>o ^ ^ then to have a balanced equation for ^o we must assume that ^,<,' ~ 

g2(i+i )_ jj^ other words asymptotic expansion near the neutral modes can only describe the soliton spectrum in 
regions of the parameter space which are close to codimension-(j + 1) bifurcation. If j =0 then only one condition 
must be satisfied and our asymptotic approach predicts presence of either two purely imaginary or two purely real 
eigenvalues, which can collide at zero. If J =1 then two conditions must be satisfied and the asymptotic approach 
predicts presence of two pairs of eigenvalues which can be real, imaginary or complex. In this situation the soliton 
becomes oscillation unstable providing that two pairs of imaginary eigenvalues collided. For each further j two new 
eigenvalues come into play. 



IV. DISCUSSION 



General formulae (|14[),(|15|) giving soliton eigenvalues with any degree of accuracy and criterion for the oscillatory 
instability ( |2l| ) are main novel results of this work. At the same time expressions for the eigenvalues near the stationary 
instability threshold, analogs of the formula A^ = Dq/Di + ..., have been earlier obtained in a number of papers. It is 
instructive now to give explicit expressions for Dj in the two simplest situations of one- and two-parameter solitons and 
to compare them with previously reported results. For the one parameter solitons: Dq = d^^Qi, Di = —{Wi\C~^Ui), 
D2 — —{Wi\C^'^Ui). Using these formulae one can show that in the case when Di ^ 0(1) the first term in Eq. ( pj| ) 
gives the same expression for A^ which was obtained in Refs. R-H], where generalised Nonlinear Shrodinger equation 
[0,yi and equations describing propagation in quadratically nonlinear media Q have been investigated. If D1D2 > 
then it can be concluded that the second term in Eq. (19) indicates saturation of the growth rate when the distance 
from the instability threshold, Dq = 0, growthes, which agrees with numerical results |qjl(|]. 



For the two-parameter solitons 
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The threshold condition Dq = has been previously found for two-parameter solitons in different physical contexts 
[QO,ni| . However derivation of an accurate expression for the soliton eigenvalues near this threshold has remained a 
controvertial problem. Indeed, comparing eigenvalues given by Eqs. (|9|), (P2|), (P3|) , ( 24 ) and eigenvalues which can 
be calculated from the ordinary differential equations for soliton parameters presented in |lQ,fl| one will discover that 

results are slightly different p^]. It has also been argued [10| that the sign of (A^j^j^ A^22 ~ AI19 -^91 )i which is the 
first term in Eq. (|2J), plays an important role in stability of two-parameter solitons. However Eqs. (p^ , (|2^) , (p3|) , ([2^) 
apparently conflict with this finding. 

Among open problems I would like to mention derivation of finite-dimensional normal forms describing dynamical 
evolution of the soliton parameters near the oscillatory instability threshold. A guideline for this work can be theory of 
G. looss p2,E3 for the normal forms of the reversible ordinary differential equations [ p3[ in vicinity of the codimension- 
2 bifurcation, wich is an equivalent of the our point Dq — Di — 0. The simplest case of the codimension-1 stationary 
instability, I?o — 0, has only one homoclinic orbit separating regions of the periodic oscillations from the spreading 
or collapse ||9|,|l3|- The vicinity of the codimension-2 point can contain the very reach dynamics, including multiple 
homoclinic orbits and stochastic regimes. 



V. SUMMARY 



General form of the asymptotic approach to stability problem of multi-parameter solitons in Hamiltonian systems 
has been developed. It has been shown that the asymptotic study of the soliton stability reduces to the calculation 



of a certain sequence of determinants, where the famous determinant of the matrix consisting from the derivatives 
of the system invariants with respect to the sohton parameters y,||,|l^ is just the first in the series. Knowledge of 
these determinants allows to calculate eigenvalues governing soliton instability with arbitrary accuracy. The most 
important consequence is that the presented approach gives first analytic criterion for the oscillatory instability of 
solitons in Hamiltonian systems. 
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FIG. 1. Soliton bifurcation diagram in the neighbourhood of the point Do = Di — for D2 > 0. Insets show 
{ReX, JmA)-plane with horizontal/vertical axes corresponding to ReX/ImX and dots marking soliton eigenvalues described 
byEq. @. 



